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ABSTRACT 



We present a measurement of the volumetric Type la supernova (SN la) rate (SNRi a ) 
as a function of redshift for the first four years of data from the Canada-France-Hawaii 
Telescope (CFHT) Supernova Legacy Survey (SNLS). This analysis includes 286 spec- 
troscopically confirmed and more than 400 additional photometrically identified SNe la 
within the redshift range 0.1 < z < 1.1. The volumetric SNRi a evolution is consistent 
with a rise to z ~ 1.0 that follows a power-law of the form (l+z) a , with a = 2.11 ± 0.28. 
This evolutionary trend in the SNLS rates is slightly shallower than that of the cosmic 
star-formation history over the same redshift range. We combine the SNLS rate mea- 
surements with those from other surveys that complement the SNLS redshift range, and 
fit various simple SN la delay-time distribution (DTD) models to the combined data. A 
simple power-law model for the DTD (i.e., oc t~^) yields values from /3 = 0.98 ± 0.05 to 
/3 = 1.15 ±0.08 depending on the parameterization of the cosmic star formation history. 
A two-component model, where SNRi a is dependent on stellar mass (M ste n ar ) and star 
formation rate (SFR) as SNRi a (z) = A x M ste \\ Sir (z) + B x SFR(z), yields the coefficients 
A = (1.9±0.1) x 10~ 14 SNeyr^M" 1 and B = (3.3±0.2) x 10~ 4 SNeyr" 1 (A^yr -1 )- 1 . 
More general two-component models also fit the data well, but single Gaussian or expo- 
nential DTDs provide significantly poorer matches. Finally, we split the SNLS sample 
into two populations by the light curve width (stretch), and show that the general be- 
havior in the rates of faster-declining SNe la (0.8 < s < 1.0) is similar, within our 
measurement errors, to that of the slower objects (1.0 < s < 1.3) out to z ~ 0.8. 

Subject headings: supernova: general — surveys 



Introduction 



Type la supernova (SN la) explosions play a critical role in regulating chemical evolution 
through the cycling of matter in galaxies. As supernovae (SNe) are the primary contributors of 
heavy elements in the universe, observed variations in their rates with redshift provide a diagnostic 
of metal enrichment over a cosmological timeline. The frequency of these events and the processes 
involved provide important constraints on theories of stellar evolution. 

SNe la are thought to originate from the thermonuclear explosion of carbon-oxygen white 
dwarfs that app roach the Chandrasekhar mass via accretio n of material from a binary companion 
(for reviews, see lHillebrandt Niemeverll2000l ; lHowellll201ll ). This process can result in a significant 



"delay t ime" between star formation and SN explosion, depending on the nature of the progenitor 
system (jMadau et al.lll998l : lGreggioll2005l ). The SN la volumetric rate (SNRi a ) evolution therefore 
represents a convolution of the cosmic star-formation history with a delay-time distribution (DTD). 
As such, measuring the global rate of SN la events as a function of redshift may be useful for 
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constraining possible DTDs and, ultimately, progenitor models - the d etailed physics of SNe la 
rema ins poorly understood, with several possible evolutionary paths (e.g. lBranch et al.lll995l ; iLivio 
20001 ) . 



One complication for rates studies is that many SN surveys at low redshifts are galaxy-targeted, 
counting discoveries in a select sample of galaxies and converting to a volumetric rate by assuming a 
galaxy luminosity function. This method can be susceptible to systematic errors if it preferentially 
samples the bright end of t he galaxy luminosity function, biasing toward SNe in more massive, or 
brighter, galaxies (see, e.g.. ISullivan et al.ll2010l ). Since many SN la properties are correlated with 
their hosts, the recovered rates may then not be representative of all types of SNe la. A second 
type of SN survey involves making repeat observations of pre-defined fields in a "rolling search", 
to find and follow SNe in specific volumes of sky over a period of time. Such surveys minimize 
the influence of host bias, but still suffer from Malmquist bias and other selection effects. It is 
reasonably straight forward — although often computationally expensive — to compensate for the 
observational biases within rolling searches. 



The advent of these wide-field rolling surveys has significantly enhanced SN statistics at cos- 
mological distances. The Supernova Legacy Survey (SNLS) in partic ular has contributed a large 
sample of Type la SNe out to redshifts of z ~ 1.05 (IGuv et alJl201(l ). Althoug 
is to assemble a sam ple of SNe la to constrain the cosmological para meters (e.g. 



1 its primary goal 



Sullivan et al 



2011 



Astier et al 
200d : 



2006; 



Bazin et al 



), the SNLS is also ideal for studies of SN rates (jNeill et al 
2009). The SNLS is a rolling high-redshift search, with repeat multi-color imaging in four tar- 
get fields over five years and as such has consistent and well-defined survey characteristics, along 
with significant follow-up spectroscopy. However, due to the selection effects (including incomplete 
spectroscopic follow-up) and other systematic errors, such as contamination and photometric red- 
shift errors, present in any SN survey, a detailed understanding of internal biases is necessary for 
accurate rate calculations. 

In the past decade, volumetric S N la rates have b e en me asured to varying degrees of accuracy 
out to redshifts of z ~ 1.6 (Fig. [1]). ICappellaro et al. ( 1999 ) compute the SN la rate in the local 
universe (z ~ 0.01) from a combined visual and photographic sample of ~ 10 4 galaxies, yet their 



by 


Li et al. (i 


2011b) 


Leaman et al. 


2011) 



using ~ 270 SNe la from the Lick Observatory Supernova Search (LOSS; 
has made signi ficant improvements in the statistics over previous studies on 



local SNe la. The rates published bv lDildav et al.l (|2010l ) include data from 516 SNe la at redshifts 



z < 0.3 from the SDSS-II Supernova Survey (SDSS-SN), with roughly half of these confirmed 
through spectroscopy. 



At intermediate redshifts, rate measurements are provided by iPain et al.l (I2002I . 38 SNe from 
the Supernova Cosm ology Project in t h e ran ge 0.25 < z < 0.85), iTonrv et al.l (|2003l . 8 SNe within 
0.3 < z < 1.2), and iRodnev fc Tonrvi feoiol . > 100 SNe from the If A Deep Survey, 23 of which 



have spectra). iNeill et al.l (|2006l ) used a spectroscopic sample of 58 SNe la from the first two years 
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Fig. 1.— Volum e tric SN la rat e s as a f unction of redshift from v arious previou s studi es , taken from 
Li et al l (hoilbh: bildav et all (|20ld ): iRodnev fc Tonrvi ifcoirt ): bahlen et all (hoosh : braur et al ' 



(2011)jjNei 



Blanc et al 



1 et a l. (2006 ). Addition a l indi vi dual rates (+) inc l ude, i n order of increasing redshift: 



(|2004l ) ; iBotticella et al.1 (]2008l ) ; iKuznetsova et al.1 ([2003) • Values are plotted as pub- 
lished, with the exception of a correction to the cosmology used in this paper. As a comp arison, the 
lines shows the evolution of vari ous model cosmic star-formation histories fromQ ( 2008, piece- wise 
fit is t he short-dashed line, the I Cole et al.l (|200ll ) form is the long-dashed line) and 
(hoosl . dot-dashed line). 



Yiiksel et al 
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of SNLS to measure a cumulative volumetric rate in the redshift range 0.2 < z < 0.6. 



SN la rates out to z ~ 1.6 are presented by IDahlen et al.l (|2004l ) using 25 SNe la (19 with 
spectra) from Hubble Space Telescope (HST) observatio ns of the Great Ob s ervato ries Origins Deep 



Survey (GOODS) fields. These data were reanalyzed bv IKuznetsova et 



identification algorithm, and the HST sample updated by IDahlen et al 



(2008) using a Bayesian 



20081 ) extending the 2004 



sa mple to 56 SNe. Groun d-based measurements from t he Subaru Deep Fie ld have also been made 
bv iPoznanski et al.l (|2007l ) using 22 SNe la, updated bv lGraur et al.l (|201ll ) with 150 events. 



The general trend of Fig. [T] reveals that the rates typically increase from z = to z = 1. There 
is a wide spread in the existing rate measurements, particularly in the range 0.4 < z < 0.8. At 
higher redshifts, data from the G OODS collabo r ation provid e some apparent evidence for a turnover 
in the SN la rates. In particular, IDahlen et al.l (|2004l . 120081 ) report a decline in SN la rates beyond 
0.8. If present, this decline mi ght point to a large r characteristic delay time between star 



z ~ 



formation and SN explosion (see also IStrolger et al.ll2004l ). However, anot her independent analysi s 
of the HST GOODS data finds ra tes that are of f set, w i th me as urements bvlKuznetsova et al.l ((2008) 
consistently lower than those of IDahlen et al.l (|2004l . |2008| ) . IKuznetsova et al.l (|2008l ) argue that 
their r esults do not distin guish between a flat or peaked rate evolution. Ground-based data in this 
range (jGraur et al.ll201ll ). while consistent with the HST-based results, show no obvious evidence 



for a decline above z ~ 1. 

In this paper we use four years of data from the SNLS sample to investigate the evolution of 
SN la rates with redshift out to z ~ 1.1. The sample presented compris es ~ 700 photomet rically 
identified SNe la from SNLS detected with the real-time analysis pipeline (jPerrett et al.ll2010l ). One 
third of these have been typed spectroscopically, and one half of the ~ 700 have a spectroscopic 
redshift (sometimes from ancillary redshift surveys in the SNLS fields) . No other data set currently 
provides such a well-observed and homogeneous sample over this range in redshift. 

Additionally, rigorous computation of the survey detection efficiencies and enhancements in 
photometric classification techniques are incorporated into the new SNLS rate measurements. 
Monte Carlo simulations of artificial SNe la with a range of intrinsic parameters are per formed 
on all of the detection images used in the SNLS real-time discovery (jPerrett et al.ll2010l ); these 
provide an exhaustive collection of recovery statistics, thereby helping to minimize the effects of 
systematic errors in the rate measurements. 

The SNLS SNe la can be used to examine the relationship between the SNRi a and redshift, 
given some model of the SN la DTD. The size of the SNL S sample also permits a division of the SNe 
la by light-curve width (in particular the "stretch"; see iPerlmutter et al.lll997l ). allowing a search 
for differences in the volumetric rate evolution expected by any changing demographic in the SN 
la population. Brighter, more slowly-declining (i.e., higher stretch) SNe la are more frequently 
found in star-forming spirals, whereas fainte r, faster- declining SNe la tend to occur in older stellar 



populations with little or no star formation (jHamuy et al 



1995 



Sullivan et al.ll2006bh . If the delay 



time for the formation of the lowest-stretch SNe la is sufficiently long (i.e., their progenitors are 
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low-mass stars ~ 10 Gyr old), these SNe la will not occur at high redshifts (jHowelll 120011 ) . The 
behavior of the high-z rates can reveal the properties of the progenitor systems. 

The organization of this paper is as follows: An overview of the rate calculation is provided 
in £J21 The SNLS data set, along with the light-curve fitting and selection cuts used to define the 
photometric sample, is introduced in £}3j SN la detection efficiencies and the rate measurements 
are presented in $3] and £}5j respectively. Several models of the SN la DTD are then fit to the rate 
evolution in £|6l and the results discussed. Finally, the stretch dependence of the rate evolution is 
investigated in £J71 We adopt a flat cosmology with (Qm ,^a)=(0.27 ,0.73) and a Hubble constant 
of H = TOkms^Mpc" 1 . 



2. The rate calculation 

The volumetric SN la rate in a redshift (z) bin z\ < z < Z2 is calculated by summing the 
inverse of the detection efficiencies, £j, for each of the N SNe la in that bin, and dividing by the 
co-moving volume (V) appropriate for that bin 

The factor (1 + Zj) corrects for time dilation (i.e., it converts to a rest-frame rate), ATj is the 
effective search duration in years, and the volume V is given by 

3 

Mpc 3 (2) 



4. e 



3 41253 



c 



z * dz' 



h j Zi ^n M (i + z>)3 + n A _ 

where O is the area of a search field in deg 2 and in this equation c is the speed of light, and Hq, 
Qm and are the cosmological parameters, and we assume a flat universe. 

Si is a recovery statistic which describes how each SN la event should be weighted relative to 
the whole population; 1 — £j gives the fraction of similar SNe la that exploded during the search 
interval but that were not detected, for example due to sampling or search inefficiencies. E{ is a 
function of the SN stretch s, a dimensionless number that expands or contracts a template light 
curve defined as s = 1 to match a given SN event, the SN color c, defined as the rest-frame B — V 
color at maximum light in the rest-frame -B-band, and the SN z. 

The Si are evaluated separately for each year and field of the survey, and are further multiplied 
by the sampling time available for finding each object (ATj) to convert to a "per year" rate. 
Typically these are 5 months for the SNLS, though this is dependent on the field and year of 
the survey. Thus, in practice, eqn. ([pj) is evaluated for each search field and year that the survey 
operates. 

This "efficiency" method is particularly suited for use with Monte-Carlo simulations of a large, 
well-controlled survey such as SNLS. Its disadvantage is that it is not straight forward to correct 
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for the likely presence of SNe that are not represented (in z/s/c parameter space) among the N in 
eqn. ([1]) (for example, very faint or very red SNe la) without resorting to assuming a luminosity 
function to give the relative fractions of SNe with different properties. In particular, we are not 
sensitive to, a nd nor do we correct for, spectroscopically pe culiar SNe la in t he SN2002cx class (e.g., 



Li et alj|2003f). and similar even ts such as SN2008ha (e.g.. lFolev et al.ll2008), super-Chandrasekha r 



Krisciunas et al. 



2011 



events (e.g.. iHowell et al.ll2006l ). and other extremely rare oddballs (e.g. 
We also exclude sub-luminous SNe la (here defined as s < 0.7, a definition that would include 
SN1991bg-like events) b ut note that these are s t udied in considerable detail for the SNLS sample 



m our companion paper. iGonzalez-Gaitan et al,l ([20 111 ). Thus, we are presenting a measurement of 
the rates of "normal", low to moderate extinction SNe la (explicitly, c < 0.6), restricting ourselves to 
the bulk of the SN la population that we can accurately model. We allow for these incompletenesses 
when comparing to other measurements of the SN la rate in which do include some of these 
classes of SNe la. 

The photometric sample begins with the set of all possible detections, to which we apply a 
series of conservative cuts to remove interlopers. The SNLS sample and the culling process are 
described next in £j3j To each resulting SN la must then be applied the corresponding e$; these 
are calculated using a detailed set of Monte Carlo simulations on the SNLS images, a procedure 
described in §3J The rate results and the measurement of their associated errors are presented 
afterwards in £|5j 



3. Defining the SNLS sample 



In this section, we describe the SNLS search and the SN la sample that we will subsequently 
use for our rate analysis. The SNLS is a rolling SN search that repeatedly t argeted four 1° x 1° 
fields (named Dl-4) in four filters (jjM r M Z M ) using the MegaCam camera (jBoulade et al.112003 ) 
on the 3.6 m Canada-France-Hawaii Telescope (CFHT). The SNLS benefited from a multi-year 
commitment of observing time as part of the CFHT Legacy Survey. Queued-service observations 
were typically spaced 3 — 4 days apart during dark/grey time, yielding ~ 5 epochs on the sky 
per lunation. Key elements of the SNLS are its consistent and well-defined survey c haracteristics. 



and the high-quality follow-up spectroscopy f rom 8m-class telescopes such as Gemini (Ho well et al 



2005 



Bronder et al 



2008; 



Walker et al. 



2011 



the ESO Very Large Telescope (VLT: iBalland et al 



2009), and Keck ([Ellis et al.ll2008l ). Due to the finite amount of follow-up observing time available, 
not all of the SN la candidates found b y SNLS were allocated for spectroscopic fo llow-up (for a de- 



scription of follow-up prioritization, see ISullivan et al. 



2006a 



Perrett et all 120101 ) . The availability 



of well-sampled light curves and color information from the SNLS nonetheless allow us to perform 
photometric identification and redshift measurements, even in the absence of spectroscopic data. 



To identify the photometric SN la sample, we begin with all variable object detections in the 
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SNLS real-time pipelin^] ( Perrett et al. 2010 ). Other articles will describe a co mplementary effor t 



to measure the rates with a re-analysis of all of the SNLS imaging data (e.g., iBazin et al.l l201ll 
We use SNLS data up to and including the fifth year of D3 observing in June, 20070. The first 
(2003) season of D3 is omitted in this analysis; this was a pre-survey phase when the completeness 
of the SN data differed significantly from the rest of the survey. The remaining detections made 
during four observing seasons for each of the four deep fields are considered in this analysis. Each 
period of observation on a given field is called a "field-season" , with 16 field-seasons in total (4 fields 
observed for 4 seasons). The coordinates of the field centers and other information are provided in 
Table ED 

We remove all candidates falling within masked regions in the deep stacks. These regions 
include areas in and around saturated bright stars or galaxies, as well as in the lower signal-to- 
noise edge regions of the dithered mosaic. The remaining unmasked areas in each field are listed 
in Table [TJ and add up to a total of 3.56 square degrees. Galaxy catalogs from these image stacks 
are used to determine the placement of test objects in the simulations described later in §H This 
cut therefore ensures that the areas being considered in the rates calculation match those used in 
the detection efficiency measurements. 

We next fit each event with a light curve fitter to determine its redshift (where no spectroscopic 
redshift is available) and photometric parameters ( §3.ip . We then remove SN la candidates with 
insufficient light curve coverage ( §3.2p . Finally, we use the light curve fits to identify and remove 
core-collapse SNe as well as other transients, such as AGN and variable stars ( §3.3f) . Each of the 
remaining SNe la will then correspond to some fraction of the true number of events having similar 
photometric properties but which were undetected by our survey. This detection efficiency will be 
determined from the Monte Carlo simulations presented in 



1 http: / /legacy.astro.utoronto.ca 

2 In June 2007, the %m filter on MegaCam was damaged during a malfunction of the filter jukebox. Candidates 
discovered after this period were observed with a new iu filter, requiring new calibrations for subsequent images, 
and were thus not included in the present study. 

Table 1. SNLS fields and survey parameters. 



Field RA (J2000) DEC (J2000) Area (sq. deg) N. 



seasons 



Dl 02:26:00.00 -04:30:00.0 0.8822 4 

D2 10:00:28.60 +02:12:21.0 0.9005 4 

D3 14:19:29.01 +52:40:41.0 0.8946 4 

D4 22:15:31.67 -17:44:05.7 0.8802 4 
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3.1. Light-curve fitting 



We fit template light curves to the SN la candidates to identify those that don't match 
typical SNe la. Flux measu rements are made on all of the final "Elixir-pr eprocessed" images 



(jMagnier Cuillandrej 12004 ) . The Canadian SNLS photometric pipeline (IPerrett et al. 



20K 



was used to measure the SN fluxes, using images processed with the accumulated flat-fields and 
fringe-maps from eac h queue run, and aligning photometrically to the tertiary standard stars of 
Regnault etafl (|200g| ). 



Two light-curve fitting tools were used to help identify the SNe la: estimate_sn (jSullivan et al 



2006 a|) for prelimin ary identification and for measuring SN la photometric redshifts, and SiFTO 
(jConlev et al.ll2008l ) for final light-curve fitting to measure the stretch and color of each candidate. 
The estimate_sn routine is not designed for exact measurement of SN la parameters, and SiFTO 
does not fit for redshift, so we require this two-step process to fully characterize the photometric 
sample of events. 

In estimate sn, t he measured fluxes in qm^m^mzm are fit using SN la spectral templates 
from lHsiao et al.l (|2007l ). The current version of the code includes the addition of priors in stretch, 
color, and Amag. These are determined from the distributions measured for the spectroscopic 
sample. The photometric redshifts (zsNphot) output from this routine are used for candidates with 
no available spectroscopic redshifts from either the SN or its host. SiFTO is an empirical light- 
curve analysis method that works by manipulat ing a spectral ener gy distribution (SED) model 
rather than operating in integrated filter space (jConlev et al.ll2008l ), SiFTO does not impose a 
color model to relate the observed filters during the light-curve fits. The implication of this is that 
SiFTO cannot easily be modified to fit for redshift, and thus requires a known input z. Output SN 
la fits are parameterized by stretch, date of maximum light, and peak flux in each filter. 

The stretch measurement provided by SiFTO is largely invariant to changes in input redshift, 
as demonstrated in Fig. [2j Even when the input redshift is off by Az = ±0.3, the output stretch 
remains within s ± 5% of its actual value. Opacit y effects in the SN ejec ta are more pronounced 
in the bluer bands, causing a more rapid decline (IKasen fc Woosleyl 120071 ) : light curves measured 
at shorter wavelengths are therefore intrinsically narrower. As a result, if SiFTO is (for example) 
given an incorrectly small input redshift, it measures the "wrong" time dilation but simultaneously 
samples further towards the blue end of the spectrum where the template is intrinsically narrower. 
The latter effect partially negates the first, resulting in the same stretch measurement regardless 
of marginal deviations in input redshift. While this is not similarly true of the derived color or 
measured fit quality, this stretch invariance is extremely useful for establishing an initial constraint 
in fitting photometric redshifts with estimate_sn. 



Spectroscopic redshifts are available for 525 (43%) of the detections remaining after the obser- 



3 CFHT-LS images processed with the Elixir pipeline are available from the Canadian Astronomy Data Centre 
(CADC): |http://www.cadc-ccd a.hia-iha.nrc-cnrc.gc.ca/cadc/ 



-0.3 -0.2 -0.1 -0.0 0.1 0.2 0.3 

Az 




s at known z 

Fig. 2. — The effects of deliberately shifting the input redshift to SiFTO. The top plot shows the 
change in output stretch for confirmed SNe la from the SNLS sample as the input redshift for the 
SiFTO fit is offset from z spcc — 0.3 to z spcc + 0.3. The gray shaded area represents the standard 
deviation of the measured points about the median As. The bottom plot shows the mean output 
stretch for each SN la as a function of its known stretch (at zero redshift offset). The error bars 
for each SN la in the lower plot represent the full range in stretch values output from SiFTO as 
the input redshift is changed. 
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vational cuts: 420 from S NLS spectroscopy and the r est from host-galaxy measurements (including 



data from DEEP/DEEP2 (jDavis et alJl2003 ). VVDS (|Le Fevre et al]l2005T ). zCOSMOS (|Lillv et al 



20071 ). and additional SNLS VLT MOS observations). The external redshifts are assigned based 
on a simple RA/DEC matching between the SNLS and the redshift catalogues, with a maximum 
allowed s eparation of 1. 5 ". For the SNLS MOS work, the host was identified following the tech- 
inques of ISullivan et al.l (l2006bl ). The known redshifts are then held fixed in the lig ht-curve fits. 



We a lso considered the use of galaxy photometric redshifts for the SNLS fields (e.g. Illbert et al 



2006). However, though these catalogs have an impressive precision, they tend to be incomplete 
and untested below a certain galaxy magnitude. SN la photometric redshifts do not suffer these 
problems. 

SN photometric redshifts (photo-z's) are calculated for the remaining objects using a multi- 
step procedure. Preliminary redshift estimates are obtained using a first round of estimate_sn 
fits without any constraints on the input parameters. The resulting fit redshifts are then used as 
input to SiFTO to measure the stretch for each object. These stretch values are then fixed in a 
subsequent round of estimate_sn fits to obtain a more robust measurement of the SN redshift - 
constraining at least one input parameter to estimate_sn improves the quality of the light-curve 
fits. Fig. [3] shows that the zsNphot are in good agreement with the spectroscopic redshifts (-z S pec) 
out to z > 0.7, with a small systematic offset above that. The median precision in zsNphot for the 
confirmed SNe la is 

MEDIAN | ; - 1 - 0.019 



{l + z t 



SpCCy 



0.031. For comparison, ISullivan et al.l (|2006al ) find |Az|/(l + z. 



''spec J 



0.031 



with cr|A^|/(i+^p C c) 

with a smaller sample and real-time data (and a previous version of the estimate_sn code). In 
§5.1| we describe how these zsNphot errors and the systematic offset are incorporated into the rates 
analysis. 



3.2. Light curve coverage cuts 

Each candidate must pass a set of light curve quality checks to be included in the photometric 
sample of SNe la for the rate calculation. Requiring that the SN light curves are well measured 
ensures that the photometric typing technique is more reliable, and that it is straight forward 
to correct for the effects of the selection cuts on the rates themselves. Therefore, candidates 
with insufficient light-curve coverage to measure accurate redshift, stretch, and color values from 
template fits are removed from the detected sample. We define observational criteria in terms of 
the phase, t, of the SN in effective days (d) relative to maximum light in the rest-frame B-b&nd, 
where 

^eff = ,° s c , (3) 
s(l + z) 

and t bs is the observer-fr ame phase of the S N. The time of maximum light is determined using the 



light curve fitter SiFTO ( jConlev et al.ll2008l ). described in the previous section 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Known z 

Fig. 3. — The comparison of SN photo- z measurements to spectroscopic redshifts for candidates in 
the SNLS sample. Filled (open) circles represent confirmed (unconfirmed) SNe la with spectroscopic 
redshifts from the SNLS, and open squares are candidates with host spectroscopic redshifts from 
the literature. 
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Each object is required to have a minimum of each of the following: 

» One observation in each of %m and tm between — 15d and +2.5d for early light-curve coverage 
and color information; 

» One observation in between — 15d and +5d for additional color information; 

» One observation in each of im and ru between — 9d and +7d for coverage near peak; 

» One observation in either iyi or ryi between +5d and +20d to constrain the later stages of 
the light curve. 



These conditions differ slightly from those used bv lNeill et al.l ((2006) in their analysis of the first 
year of SNLS data. Note that no cuts are made on the signal-to-noise (S/N) on a particular epoch; 
that is, a detection of a candidate on each of the observation epochs is not a requirement. We also 
neglect the redshift offset seen in Fig. [3] in calculating the above rest-frame epochs. We estimate 
that this would shift the effective epochs by only one day in the worst case (+20d; a z = 1 SN), 
and in most cases would be far smaller than this. 

Table [2] provides the numbers of candidates that survive each of the applied cuts. In total, 1210 
SNLS detections pass the light curve coverage cuts, 305 of which are spectroscopically confirmed 
SNe la. (For consistency, these same objective requirements are also applied to the artificial SNe la 
used in the Monte Carlo simulations (Q, thereby directly incorporating the effects of this cut into 
the detection efficiencies.) With these objective criteria satisfied, we can then use light-curve fitting 
to define a photometric SN la sample. 



3.3. Removing non SNe la 

A set of xt goodness-of-fit cuts (here, xt is the x 2 per degree of freedom, v) are applied to all of 
the SN la light-curve fits from estimate_sn to help eliminate non-las from the current sample (see 



also ISullivan et a l. 2006a). An overall xt cu t along with individual %m and tm filter xt constraints 



are applied separately for cases both with and without known redshifts. Light-curve fit quality 

Table 2. Number of candidates after selection cuts 



Cut 


All candidates 


Confirmed SNe la 


Masking cut 


1538 


325 


Observational cuts 


1210 


305 


Fit quality and s cuts 


691 


286 
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limits for the sample with known input redshifts are set to Xu < 9, x\ < 9> an d Xr < 18 (here the 
is omitted for clarity). Those with fit redshifts are given stricter limits of Xu < 6, x| < 6, and 
< 12. The tighter xt limits for candidates without known redshifts are necessary since core- 
collapse SNe — specifically SNe Ib/c — can sometimes achieve better fits to SN la templates when 
z is permitted to float from the true value. The limits are determined empirically by maximizing 
the fraction of SNe la remaining in the sample, while also maximizing the number of known non- 
las that are removed. Note that estimate_sn, unlike SiFTO, enforces a color relation between the 
fluxes in different filters, which leads to larger xt than in SiFTO fits. 

In the case of fixed [floating] redshifts input to the fit, > 95% [> 96%] of spectroscopically 
identified SNe la survive the xt cuts (we correct for this slight inefficiency when calculating our 
final rate numbers), while 0% [13%] of confirmed non-las remain in the sample. A final round 
of SiFTO fits is then performed to determine the output values of stretch and color. The input 
redshifts to SiFTO are set to zsNphot wherever no z spec values are available. 

One final light-curve fitting cut is then applied on the sample, requiring that the output SiFTO 
template fits have XsiFTO ^ 4. This step removes all but one of the remaining confirmed non-Ias3 
when all redshifts are allowed to float, while at the same time maximizing the number of confirmed 
SNe la passing the cut. No known contaminants remain when all available z spec values are fixed in 
the fits. 



3.4. The photometric SN la sample 

The final photometric SN la (phot-la) sample is restricted to 0.1 < z < 1.1. Above this 
redshift, the rates are found to be too uncertain to include in subsequent analyses. This is a result 
of low S/N, poor detection efficiency, 100% spectroscopic incompleteness, and the potential for 
increased contamination from non-las. Only candidates having stretch values within 0.7 < s < 1.3 
are considered in the present study. This range is characteristic of the SNLS spectroscopic sample 
- shown by the red histogram in the central plot of Fig. U] - but excludes sub-luminous events 
such as SN1 991bg. These sub-luminous, lo w-stretch SNe la in the SNLS sample have been studied 



in detail by iGonzalez-Gaitan et al.l (|201ll ) - our stretch limit removes 22 such objects from our 
sample. Extremely red (c > 0.6), and presumably highly extincted, candidates are also removed. 
This cut eliminates only one event: SNLS-04D2fm, a faint SN of unknown type at z spcc = 0.424. 

The final redshift, stretch, and color distributions resulting from the various cuts are shown in 
Fig. SJ The phot-la sample consists of 691 objects, 371 of which have known redshifts. A total of 
286 objects in this sample have been spectroscopically confirmed as Type la SNe (Table [2]). The 
redshift histogram reveals that the incompleteness of the spectroscopic sample (in red) begins to 



4 The identification of SNLS-06D4cb is inconclusive, although it has a spectrum that is a poor match to a SN la. 
The SN photometric redshift for this object is Zfu = 0.64 but the host has a spectroscopic redshift of z = 0.4397. 
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Redshift (z) 
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Known z [I 
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Color (c) 
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Stretch (s) 



1.4 



Fig. 4. — Distributions in redshift (upper left), stretch (upper right), and color (lower left) for the 
SNLS SNe la. The gray histogram represents the final photometric SN la sample and the blue 
histogram shows the fraction of the sample with known redshifts. The spectroscopically confirmed 
SNe la are shown as the red distribution in each plot. Sample incompleteness causes the decline in 
the observed population at z > 1.0. 
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increase beyond z ~ 0.5, where the rise in the required exposure time makes taking spectra of 
every candidate too expensive. The effects of incompleteness in the observed SNLS sample become 
severe above z > 1.0. The full phot-la sample has median stretch and color values of s = 1.00 and 
c = —0.04, respectively. The color distribution peaks at a slightly redder value than the estimated 
typical color of a SNe la of Cf ~ —0.06, based on the distribution observed for the spectroscopic 
SNLS sample. 



4. Detection efficiencies 



With the final SN la sample in hand, we now need to estimate the weight that each of these 
events contributes in our final rates calculation. These "detection efficiencies" depend on many 
observational factors and will obviously vary with SN la characteristics. For example, at higher- 
redshift, the higher stretch SNe la are more likely to be recovered not only because they are brighter, 
but also because they spend a longer amount of time near maximum light, and are therefore more 
likely to pass the culls of In a rolling search like SNLS, such effects can be directly accounted 
for by measuring recovery statistics for a range of simulated input SN la properties using the 
actual images (and their epochs) observed. This is a brute-force approach, but is a practical way to 
accurately model a survey such as SNLS, helping to control potential systematic errors by avoiding 
assumptions about image quality limitations and data coverage that may bias the rate calculation. 
Uncertainties on search time and detection area are avoided since the actual values are well defined. 



4.1. Monte Carlo simulations 



An exhaustive set of Monte Carlo simulations were performed for each field-season to determine 
the recovery frac tion as a funct i on of redshift, stretch, and color. Full details about these simulations 
are presented in lPerrett et al.l ( 201C ). 

A total of 2.5 x 10 6 artificial SNe la with a flat redshift distribution were added to galaxies 
present in the SNLS fields. Each host galaxy was chosen to have a photometric redshift within 
0.02 of the artificial SN redshift, with the probabil ity of selecting a parti cular gala xy we ighted by 
the "A+B" SN rate model with coefficients from I Sullivan et al.1 (|2006bl . hereafter |S06| ) . Within 
their host galaxies, the artificial SNe were assigned galactocentric positions drawn from the two- 
dimensional Gaussian distribution about the host centroid returned by SExtractor, i.e., the artificial 
SNe are placed with a probability that follows the light of the host galaxy. 

The simulated objects were assigned random values of stretch from a uniform distribution in 
the range 0.5 < s < 1 . 3, wit h colors calculated from the stretch-color relationships presented in 
Gonzalez-Gaitan et al.1 (|201ll ) (the use of a uniform distribution in stretch ensures that the pa- 
rameter space of SN la events is equally sampled). Peak apparent rest-frame B magnitudes (tub) 
at each selected redshift were calculated for our cosmology and a SN la absolute magnitude, and 
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adjusted for the color-luminosity and stretch-luminosity relations. We use an em pirical piece-wise 
stret c h-luminosity relationship wit h different slopes above and below s = 0.8 (e.g. iGarnavich et al 



2004 



( Conlev et al 



2011 



Sullivan et al 



Gonza lcz-Gait an et al .1 12011 ) , and SN la photometric parameters from the SNLS3 analysis 



20 111 ). These peak apparent magnitudes were then further ad- 



justed by an amount Amag according to the observed intrinsic dispersion (<7i nt ) in SN la magnitudes 
following s and c corrections. Here, <7i n t parameterizes a Gaussian distribution from which a Amag 
can be assigned for each artificial event. 

The SN color-luminosity relation includes both effects intrinsic to the SN, and extrinsic effects 
such as dust. We use coefficients consistent with the SNLS3 analysis, which favor a slope between 
tub and c of < 4.1, the valu e expected based on Milky Way dust. As there is no evidence that this 



slope evolves with redshift (jConley et al. 



2011 



), we keep it fixed for all the artificial SNe. For the 
detection efficiency grids, our c values range up to 0.6, corresponding to a SN that is ~ 1.8 mag 
fainter in £?-band than a normal SN la. 

Each artificial SN was assigned a random date of peak magnitude. For the field-season under 
study, this ranged from 20 observer-frame days before the first observation, to 10 days after the 
last observation. This ensures that the artificial events sample the entire phase range allowed by 
the culls in § 13.21 at all redshifts. The light curve of each event in %m was then calculated using the 
^-correction appropriate for each epoch of observation, and each artificial object was added a t the 
appropriate magnitude into every %m image. The real-time search pipeline (IPerrett et al.ll2010i ). the 
same one that was used to discover the real SNe, was run on each epoch of data to determine the 
overall recovery fraction as a function of the various SN la parameters. The variation in candidate 
recovery over magnitude, redshift, stretch, and color are shown by the solid lines in Fig. [5j The 
50% detection incompleteness limit lies at im = 24.3 mag in the AB system. 

As expected, SNe la that are high stretch, blue, or at lower redshift are all generally easier to 
recover. Note that at lower redshifts, the faster (less time-dilated) nature of the SN la light curves 
means that the observational criteria of §3.21 are slightly more likely to remove events (as there are 
fewer opportunities to observe a faster SN), hence the observed decrease in the recovered fraction 
towards lower-redshifts. That is, a low-z SN that peaks during bright time is less likely to be 
recovered than a higher- z SN peaking at the same epoch, even if they had the same observed peak 
magnitude. This is also partially reflected in the fraction recovered as a function of magnitude, 
with a curvature in t he recovered fractio n towards brighter magnitudes. The recovery results are 
discussed in detail in lPerrett et all (|2010l ). 



A grid of detection efficiencies was constructed independently for each field-season using the 
recovery statistics in bins of measured redshift (Az = 0.1), stretch (As = 0.1), and color (Ac = 0.2). 
These bin sizes were found to provide adequate resolution in each parameter. We investigated the 
use of a higher resolution in stretch and color, and found no significant on impact our results. 
Every observed phot-la in the SNLS sample is thereby assigned a detection efficiency by linearly 
interpolating in z/s/c space that corresponds to the field-season during which it was detected, 
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Fig. 5. — Recovery fraction as a function of i\j (AB) magnitude (upper left), redshift (upper right), 
stretch (lower left), and color (lower right) for all field-seasons combined. The solid lines represent 
the fraction of objects found, and the dashed lines include the additional observational constraints 
as described in §3.21 These plots include only the artificial SNe la from the simulations that lie 
within the parameter space typical of the observed SNLS sample. 
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along with its other measured parameters: e(field, z, s, c). These detection efficiencies are plotted 
in Fig. [6] prior to any adjustments for sampling time and the availability of observations. Redder, 
lower-stretch SNe la tend to have smaller detection efficiencies, as shown by the open circles in 
Fig. [6l For clarity, detection efficiency errors are not shown in Fig. [6l 

Statistical uncertainties on e(field, z, s, c) for well-sampled data are governed by the number of 
Monte Carlo simulations performed, and are small in comparison to the systematic error resulting 
from assumptions made about the underlying intrinsic SN la magnitude dispersion (the Amag 
distribution, parameterized by <7i n t). To estimate these latter errors, the detection efficiencies are 
recalculated using a range of <7i n t values from 0.12 — 0.15. Bins with e = 1 will effectively have zero 
uncertainty, since the likelihood of recovery will not depend on the details of the population dis- 
tribution; by contrast, "low-efficiency" bins are more seriously affected. These detection efficiency 
"errors" are included into the overall rate uncertainties in §5.11 



4.2. Sampling time 

To remain consistent in the selection criteria used for both the observed SNLS sample and 
the fake objects, we also apply the same observational cuts described in §3.21 to the artificial 
SNe la. Using the peak date of each simulated light curve, we determine whether the minimum 
observing requirements are met in each filter by comparing with the SNLS image logs. This directly 
incorporates the observational cuts into the detection efficiency calculations, while factoring in losses 
due to adverse weather and the gaps between epochs. The recovery fractions that include these 
observational requirements are shown by the lower dashed lines in Fig. [5j 

Each candidate's detection efficiency is multiplied by a factor to account for its corresponding 
sampling time window for detection, yielding a "time corrected" rest-frame efficiency et- 

1 AT /j\ 

£ T = e—— (4) 

(1 + z) yr 

The sampling period AT (in years) for a given field-season is 

AT = 365 1 ^ [max(MJD) - min(MJD) + 30] , (5) 

where MJD is the modified Julian date of the available detection images. The extra 30 days account 
for the range in peak dates allowed for the artificial SN la light curves, from 20 days prior to the 
first observation in a given field-season to 10 days past the final epoch. 

The resulting time-corrected rest-frame detection efficiencies for the phot-la sample are plotted 
as a function of redshift in Fig. UJ Since each field is observable for at most 4 — 6 months of the 
year, et peaks at ~ 0.4 even for bright, nearby objects. 

Figs. [6] and [7J show that there is a drop-off in the efficiencies above z = 0.9, in particular for 
the redder c bins, making it more difficult to calculate accurate rates at these redshifts due to 
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Redshift 

Fig. 6. — Detection efficiencies (e, in eqn.Q]) measured for each candidate in the photometric SN la 
sample and plotted against redshift. SNe la that are redder than the adopted fiducial color of 
Cf ~ —0.06 are shown as red circles, while bluer objects are shown as blue squares. Open symbols 
represent SNe la with stretches smaller than the median value of the sample (s < 1). These 
efficiencies have not been corrected for changes in the sampling time between the different fields 
observed. 
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Redshift 

Fig. 7. — Time-corrected rest-frame efficiencies for the SNLS phot-la sample plotted against red- 
shift. The efficiencies shown here are < 1 even at low-z since they have been adjusted for field 
observability. The dashed line shows that a (1 + z)~ l slope matches the general trend of the data 
out to z ~ 1, where the detection efficiencies begin to drop off more quickly. SNe la with c > Cf 
are shown as red circles and bluer ones as blue squares. Lower-stretch (s < 1) events are displayed 
as open symbols. There are no significant differences in the median values of £t as a function of 
redshift for high- and low-stretch objects. 
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color-stretch bins which are not sampled. At z > 1.1, it is not possible to measure SN la rates 
using this method due to poor survey sensitivity and inadequate statistical sampling of spectral 
templates. Therefore, we restrict our volumetric rate calculations to the range 0.1 < z < 1.1. 

5. SN la rates 

Volumetric SN la rates are calculated from eqn. Q] by summing the observed SNe la weighted 
by the inverse of their time-corrected rest-frame efficiencies. The total sampling volumes for the 
deep fields in redshift bins of Az = 0.1 (eqn. [2D are provided in Column 2 of Table EJ Columns 3 
and 4 show the numbers of observed candidates in each bin for the entire sample (iV b s ) and for the 
spectroscopically confirmed SNe la (iV S p ec _i a ) in each redshift bin. The "raw" measured rates (r meas ) 
with their weighted statistical errors are given in Column 5, in units of xlO -4 SNeyr -1 Mpc -3 ; see 
later sections for the meaning of the remaining columns. 

Contamination by non-las that survive the culling criteria is estimated to contribute under 2% 
to the total measured rates to z ~ 1. The contribution is found to be negligible up to z ~ 0.5, 
at which point it increases to around 4% at z ~ 1. This is determined by summing 1/er for 
the known non-las, and dividing by the corresponding value for objects in each redshift bin with 
available spectroscopy. The £t values used here are based on the results obtained when allowing 
redshift to vary in the fits, not when holding z fixed at the spectroscopic redshifts. 

We now correct the raw measured rates for potential systematic offsets in the photometric 
redshifts and other parameters. This is done using the technique described next in §5-H which also 
computes a combined statistical and systematic uncertainty on the final rates. The SN la rates are 
potentially also sensitive to the inclusion of very low detection efficiency candidates at z > 0.9, and 
we must consider the effects of undetected SNe la in z/s/c bins with very poor detection recovery 
rates. These low-efficiency issues are discussed later in £ j5,2i 

5.1. Error analysis 

In addition to the simple "root-iV" statistical errors, a number of additional uncertainties also 
affect our measured rates. These can include errors in the measured SN (photometric) redshift, 
stretch, and color, which together determine the detection efficiency (and hence weight) assigned 
to each event. 

As shown in Fig. [3j there is a redshift uncertainty in each measure of ^sNphot • The discrepancy 
between z spec and ^sNphot for the confirmed SN la sample is presented in Fig. [8] for two bins in 
stretch. There is a small offset above z = 0.7, increasing to Az = z spec — zsNphot 

fa 0.05 {a = 0.08) 

at z > 1 in the 0.7 < s < 1.0 sample, and 0.07 (cr = 0.06) in the 1.0 < s < 1.3 sample. On 
average, the zsNphot measurements are underestimated, with an increasing offset to higher redshift. 
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An offset is expected based on a Malmquist bias, such that bri ghter objects are more likely 
to have a spectroscopic type at a fixed redshift (jPerrett et al.ll2010l ). However, we estimate this 
effect to be smaller: the solid lines in Fig. [8] show this predicted offset as a function of redshift. 
This is calculated by estimating the rest-frame U-band apparent magnitude with z for the adopted 
cosmology, and app lying the Amag offsets contributed by spectroscopic selection as measured in 
Perrett et all (bold ). 



To study these various uncertainties, and to handle this redshift migration effectively, we 
perform a set of Monte Carlo simulations on the measured rates. We begin with the basic rate 
measurements, r meas (z), from Tabled and calculate how many measured SNe la that rate represents 
in each redshift bin by multiplying by the volume in that bin: N meas (z). Many realizations (5000) 
are performed by drawing iV meas objects from typical SNLS-like distributions of artificial SNe la. 
These are the same artificial objects as used in the detection efficiency calculations (Si), although 
the stretch, color, and Amag distribut ions are matched to those of the spectroscopically confirmed 
SN la sample (see iPerrett et al.ll2010l ). The distributions for input to the error calculations are 
shown in Fig. [9j Only a fraction of the objects in each redshift bin have a spectroscopic redshift, 
with the remainder having a SN la photometric redshift and accompanying uncertainty (Fig. [8]). 
This "spectroscopic fraction", F spec (z) is calculated in each redshift bin from Fig. UJ 



The procedure for each realization is then as follows: 



1- ^mcasl 2 ) is randomized according to the Poisson distribution, using the Poisson error based 
on N b s (z) but scaled to N meas (z), to give N ian ^(z) simulated objects. 

2. Each simulated object i in each redshift bin is assigned a random redshift appropriate for that 
bin (zi). Within each bin, the probability follows a scaled number-density profile according 
to the expected increase in volume with redshift. 

3. The Zi are then randomly matched to a SN la from the artificial distribution with the same 
redshift (Fig. [9|), and that event's stretch (sj) and color (cj) assigned to the simulated event. 

4. Using the fraction of spectroscopic redshifts in each bin F spec (z) (Fig. |4]), we assign this 
fraction of the -/V ran d objects to correspond to a spectroscopic redshift measurement. The 
remaining redshifts are assumed to come from a photometric fit, and are shifted and random- 
ized using the median offsets and standard deviations shown in Fig. [8] to give z[. Any event 
with a z spec is not adjusted. 

5. Correlated stretch and color errors are then incorporated for all objects, using typical co- 
variances produced by SiFTO for the SNLS sample, and the s, and c, randomized to and 

6. z[, s[, and c- are used to match each simulated object to a detection efficiency. Efficiency 
errors are included by applying a random shift drawn from an two-sided Gaussian representing 
the asymmetric uncertainties on each value ( §4.11) . 
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Table 3. Volumetric rates from the SNLS sample 



z bin Survey Volume V 

(10 4 Mpc 3 ) 



Nn 



N. 



spec — la 



r 

(xlO 



meas ' mcas 

■^SNeyr^Mpc" 3 ) 



(xlO 



^SNeyr-^pc" 3 ) 



0.10 


-0.20 


17.3 


4 


3 


0.21 ±0.11 




0.16 


0.14 


0.20 


-0.30 


42.8 


16 


16 


0.30 ±0.08 




0.26 


0.28 


0.30 


-0.40 


75.7 


31 


24 


0.35 ±0.07 




0.35 


0.36 


0.40 


-0.50 


112.7 


42 


29 


0.36 ±0.06 




0.45 


0.36 


0.50 


-0.60 


151.5 


72 


47 


0.48 ± 0.06 




0.55 


0.48 


0.60 


- 0.70 


190.1 


91 


36 


0.55 ± 0.06 


0.57 ±0.06 


0.65 


0.48 


0.70 


- 0.80 


227.2 


110 


56 


0.59 ± 0.06 


0.57 ±0.06 


0.75 


0.58 


0.80 


- 0.90 


262.1 


128 


44 


0.64 ± 0.06 


0.65 ± 0.06 


0.85 


0.57 


0.90 


- 1.00 


294.1 


141 


25 


1.20 ± 0.17 


0.99 ± 0.29 


0.95 


0.77 


1.00 


- 1.10 b 


323.0 


50 


6 


0.93 ± 0.25 


0.51 ±0.26 


1.05 


0.74 



,,+0 

-0 
2 + 

<t8: 

q+O. 

s+o 
a+o 

J+O 

7+0 
-0 
7+0 
-0 
1+0 



09+0.06 
09-0.12 
07+0.06 

07- 0.07 
06+0.05 
06-0.06 
06+0.04 
06-0.05 
06+0.04 
06-0.05 
05+0.04 

05- 0.06 
06+0.05 

06- 0.07 
05+0.06 
05-0.07 
08+0.10 

08- 0.12 
12+0.10 
12-0.13 



a The first error listed is statistical, and the second systematic. 

b Bins at z > 1.0 are not included in the rates analysis; see Section \5. 21 
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Fig. 8. — The redshift offset, Az = z spec — zsNphot > as a function of SN or host spectroscopic redshift 
for the phot-la sample. The offsets are calculated separately in two stretch bins: 0.7 < s < 1.0 
(upper panel) and 1.0 < s < 1.3 (lower panel). The median z offset in sliding bins of width Az = 0.2 
are shown by the solid points, with error bars representing the standard deviation in each bin. The 
offsets increase from approximately zero at z = 0.7 to Az ~ 0.05 — 0.07 at z > 1. The solid lines 
represent the expected offset due purely to sample selection bias in each stretch range. 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.6 0.8 1.0 1.2 1.4 
Redshift Stretch 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 
Color Amag 

Fig. 9. — Resampled distributions showing the properties of the artificial SNe la used as input to 
the rate error simulations. Amag refers to the scatter in SN la rest-frame 5-band peak magnitudes, 
and has a dispersion of fTi Ilt = 0.14. 
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Fig. 10. — Histograms showing the input (black line) and output (gray filled) parameter distri- 
butions from the rate error simulations as a function of redshift (left), stretch (center), and color 
(right). 
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7. Random numbers between zero and one are generated to evaluate whether each simulated 
object gets "found": if the selected number is lower than the detection efficiency associated 
with the simulated object, that event is added to the rate calculated for that iteration. 

8. The rate for that iteration is then calculated using the appropriate detection efficiencies from 
step 6. 

The final volumetric rates are presented in Columns 7 of Table [3l These are calculated as the 
mean of the 5000 simulated rates in each redshift bin. We also calculate the standard deviation 
in each bin, subtract from this in quadrature the statistical Poisson uncertainty based on N ^ s (z), 
with the remainder our estimate of the systematic uncertainty in each redshift bin. 

The effects of the simulations described above on the input redshift, stretch, and color distri- 
butions are shown in Fig. [10l The redshift histogram shows that the offsets applied in step 4 of the 
simulations produce a net increase in redshift to compensate for the small bias in the photometric 
fitting, with the effect increasing towards higher z. This causes a flattening of the output rates 
calculated by the simulations as compared with the measured (and uncorrected) rates (Table [3]). 
In addition to the offset towards higher redshifts at z > 0.7, there is also a very small spread in the 
stretch and color distributions (Fig. fTUl) . 

5.2. Low-efficiency candidates 

The detection efficiencies in some of the reddest c bins begin to rapidly decrease at z > 0.9 
(Fig. [7]). As the contribution to the rate from each observed SN la goes as 1/et, the measured 
volumetric rates are particularly sensitive to any objects with very low detection efficiencies (er)- 
There is also the potential for a complete omission of SNe la in some z/s/c bins. For example, in 
z/s/c bins with less than 10% detection efficiency, on average at least 10 SNe la must be observable 
in a given field-season for just one to be detected. If that one SN is not detected by the real-time 
pipeline, the 10 SNe la that it truly represents in that bin will never be counted in the final rates 
tally (and the rate measurement will be biased). 

To examine the sensitivity of our rates on the "low-e^ regions" of z/s/c parameter space, we 
use the z < 0.6 detection-efficiency-corrected SNLS sample as a model for the true (s, c) SN la 
distribution at z > 0.6. This population is assumed observationally complete (Fig. [5]) and is taken 
to be representative of the underlying sample of SNe la in the universqj. The two-dimensional (s, 
c) distribution at z < 0.6 is fit to the five z > 0.6 bins, and the best-fit scaling determined. The 
total rates r' meas (z) are then calculated from those scaled numbers (tabulated in Table [3]) 



Of course, this relies on the (possibly incorrect) assumption of no evolution in intrinsic stretch or color as a 
function of redshift (e.g.. iHowell et al.1 2007). 
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These tests indicate that, while the results remain consistent within their errors up to z = 0.95 
(0.99 ± 0.29 x lO^SNeyr^Mpc" 3 compared with r v = 1.20 ± 0.17 x 10~ 4 SNeyr^ 1 Mpc~ 3 ), 
there is a significant amount of uncertainty in the SNLS rates at higher redshifts due to sample 
incompleteness. At z = 1.05, the scaled rate is 0.51 ± 0.26 x 10 -4 SNeyr -1 Mpc~ 3 , whereas our 
the calculated value is ry = 0.93 ± 0.25 x 10 -4 SNeyr -1 Mpc -3 . This finding is consistent with the 
results shown in Figs. |4] and El the phot-la sample numbers drop significantly at z > 1.0, and those 
that are found in the sample can have very low ey. For these reasons, we limit the formal analysis 
of the SNLS rates to r v (z < 1.0). 



6. Delay-time distributions 



Having measured volumetric SN la rates and associated errors, we now compare our mea- 
surements with those of other studies. We also examine the SN la rate evolution as a function 
of redshift, and compare with predictions based on various simple delay-time distribution models 
from the literature. 

For comparison of the SNLS rate measurements to various SN la models, additional data in 
redshift ranges not sampled by SNLS is requir ed. In the rest o f this section, we will make use of 
an extended SN la rate sample comprising the iLi et al.l (|2011bl ) LOSS measuremen t at z ~ 0, the 
Dildav et al.l (|2010l ) sample from SDSS-SN at z ~ 0.2, and the recent iGraur et al.l (|201ll ) Subaru 
Deep Field (SDF) sample at higher redshifts, together with our SNLS results. Clearly other samples 
could have been chosen - however, these three are the largest SN la samples in their respective 
redshift ranges, and have the greatest statistical power. In the case of the SDSS and SDF samples, 
they are also built on rolling SN searches similar to SNLS. 

We make some s mall corre c tions t o these published rates in order to ensure a fair comparison 
across samples. The ILi et al.l (|2011bl ) sample includes all sub-classes of SNe la, including the 
peculiar events in the SN2002cx-like class and sub-luminous events in t he SN1991bg-lik e class. 
SN2002cx-like events make up 5% of the LOSS volume-limited SN la sample (|Li et al.ll2011al ). These 
are not present (or accounted for) in the SNLS sample, a nd are ex c luded f rom the SDSS analysis 
(jDilday et al.ll2008l ). so we therefore exclude these from the lLi et al.1 (j2011bl ) sample, reducing their 
published rate value by 5%. 



Both the lLi et al.l (|2011bl ) and Dilda v et al.l (20101 ) samples include SNe la in the sub-luminous 
SN1991bg category (s ee also iDildav et al.l 120081). which we exclude here in the SNLS analysis 
(these are studied in iGonzalcz- Gaitan et al.l 120111). While we could correct our own rates for 



this population using the iGonzalez-Gaitan et al.l (|201ll ) results, it is unclear how to treat the SDF 



sample in the same way (do SN1991bg-like events even occur at z > 1?) and the SNLS sub-luminous 
measurement is quite noisy. Instead we use the very well-measured fraction of SN1991bg-like SNe 
in the volume-limited LOSS sample (15%), and reduce both the L OSS and SPSS p ublished rates 
by this amount. This 15% is based on the classifications given in ILi et al.l (j2011bl ). We confirm 
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that this is app ropriate for our stretch selection (i.e., we req uire s > 



Li et al.l (|2011bl ) light curves with SiFTO. 16% of the available iLi et ajj (2011b) sample has a fitted 



stretch < 0.7, consistent with the 15% reported as 91bg-like bv lLi et al. 



7) by fitting the available 



(|2011bT ). 



6.1. Comparison with published rates 



Fig. [TT] shows the SNLS volumetric SN la rates f or comparis o n wit h recent published results. 
The SNLS volumetric rate at z ~ 0.5 published by iNeill et al.l faood ) is SNRi a «z) = 0.47) = 
[0.42l°; 43 (syst) ± 0.06(stat)] xl0~ 4 SNeyr" 1 Mpc" 3 , consistent with our binned rates in the same 
redshift range. Note that the effects of non-la contamination (£}5]) have not been incorpo r ated into 
the SNLS errors shown in Fig. 111! Ou r results are also consistent with iDildav et al.l (|2010l ) at 
< 0.3, and with iRodnev fe Tonrvi (|2010l ) at higher redshifts, although those latter measurements 
have significantly greater uncertainties. 

The SNLS rates show a rise out to z ~ 1, with no evidence of a rollover at z ~ 0.5. We can 
parameterize the SNLS rate evolution as a simple power-law: 

SNR Ia (*) = r (l + z) a , (6) 

with the best-fit shown as the solid line in Fig. [[2j We find a = 2.11 ± 0.28 and r = (0.17 ± 
0.03) x lO^SNeyr^Mpc" 3 (all errors in this section are statistical only). This evolu tion is 
shallower than a typical fit to the cosmic star-formation history, with a ~ 3.3 (e.g. [Ei 2008 ). 
though the constra ining power of the SNLS data alone at z < 0.3 is not great. By comparison 



Dildav et all (hold ) find r = (0.23 ± 0.01) x 10" 4 SNeyr" 1 Mpc" 3 and a = 2.04±°J° using the 



lower-redshift SDSS-SN data, completely consistent with our results. Including all the external 
data gives r = (0.21 ±0.01) x 10~ 4 SNeyr" 1 Mpc~ 3 and a = 1.70 ±0.12. 



6.2. Comparison with delay-time distribution models 



We now compare our rate evolution with simple parameterizations of the delay-time distribu- 
tion (DTD) from the literature relating the cosmic star-formation history (SFH) to SN la rates. 
The DTD, ^(t), gives the SN la rate as a function of time for a simple stellar population (SSP), 
i.e., following a 5- function burst of star formation. The SN la rate (SNRi a ) at time t is then 



SNRia(i) 



SFR(t - r)^(r)dr. 



(7) 



where SFR(i) is the star-formation rate as a function of time. Thus, different functional forms of 
the DTD can be tested against observations of volumetric rates if the SFR(t), or the cosmic SFH , 



is known (e.g.. Madau et al.lll998l : IStrolger et al.ll2004l : lOda et al.ll2008l : iHoriuchi fc Beacomll2010l : 



Graur et al.ll201ll ). An implicit assumption in this test is that the DTD is invariant with redshift 
and environment. 
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Fig. 11. — The SNLS volumetric SN la rates in the context of the data in Fig. [TJ The filled circles 
represent the SNLS rates from the current analysis. The rate at z = 1.05 (with the dashed error 
bar) represents the redshift bin in which incompleteness and poor spectroscopic sampling make 
measureme nts untrustworth y. S NLS rates abo v e z = 1.0 are not included in subsequent fits. The 
samples of iLi et al.1 (|2011bl ) and iDilday et al.1 (j2010l ) have been scaled downwards to reflect the 
exclusion of sub-luminous and SN2002cx-like SNe la from the SNLS sample (see §[6]). Over-plotted 
are the various SFHs we use in our anal ysis in § [6] as fit to the SNLS data only. The short-dashed 
l ine sh ows the piece-wise SFH from Li ( 2008), the long- d ashed line the Cole et al. form from Q 
(boosh. the dot-dash ed line the SFH from lYiiksel et al.1 l|2008l ). and the dotted line the SFH of 
Wilkins et all rt2008h . 
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Fig. 12. — SNLS rates as a function of redshift, showing a power-law fit to the data (solid line): 
SNRia(z) = r (l + z) a , where a = 2.11 ± 0.28 and r = (0.17 ± 0.03) x 10~ 4 SNeyr" 1 Mpc~ 3 . 
The reduced y 2 goodness-of-fit st atistic is xi = 0.64. For comparison, the dashed line shows the 
Cole et al. ( 2001 ) form of the Q ( 2008 ) star formation history profile, which has a = 3.3 out to 



z ~ 1. 
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Table 4. Various DTD model fits to the volumetric SN la rate, SNRi a . 









Ij (2008) piece-wise 


SFH 


Li (2008) "Cole et al." 


SFH 


Yiiksel et al. (2008) SFH 


Data 


Model 


V 


Fit parameters 


X 1/ 


Fit parameters 


X.i/ 


Fit parameters 










r (Gyr) 






t (Gyr) 






r (Gyr) 






Extended 3 


Gaussian 


18 


3.4 




2.62 


3.4 




1.60 


3.4 




3.29 


Extended 


Gaussian 


17 


3.1 ±0.3 




2.64 


2.5 ± 0.7 




1.60 


2.9 ±0.3 




3.00 


Ext.+D08 


Gaussian 


19 


3.4 




2.72 


3.4 




1.68 


3.4 




3.43 








P 






P 






P 






SNLS 


Power law 


8 


1 




0.76 


1 




1.36 


1 




0.81 


Extended 


Power law 


18 


1 




0.72 


1 




1.06 


1 




0.78 


Extended 


Power law 


17 0.98 ±0.05 




0.76 


1.15 ± 0.08 




0.92 


0.98 ±0.05 




0.81 








r (Gyr) 






r (Gyr) 






r (Gyr) 






SNLS 


Exponential 


7 


1.5 ±0.4 




0.85 


0.2 ± 2.9 




0.58 


1.4 ±0.4 




0.91 


Extended 


Exponential 17 


2.6 ±0.3 




1.33 


2.1 ± 0.3 




1.15 


2.5 ±0.3 




1.44 


P P P 


SNLS 


PHS 


8 


0.5 




4.14 


0.5 




4.81 


0.5 




4.31 


Extended 


PHS 


18 


0.5 




4.08 


0.5 




4.77 


0.5 




4.19 








A a 


B h 




.4 


B 




.4 


B 




SNLS 


A + B 


7 


1.6 ±0.5 


3.4 ±0.3 


0.74 


0.3 ±0.6 


5.2 ±0.5 


0.59 


1.7 ±0.5 


3.2 ± 0.3 


0.78 


Extended 


A + B 


17 


1.9±0.1 


3.3 ±0.2 


0.60 


1.5 ± 0.2 


4.3 ±0.3 


0.77 


2.0 ±0.1 


3.1 ± 0.2 


0.63 


SNLS 


A = 


8 





4.4 ±0.3 


1.81 





5.4 ±0.2 


0.54 





4.2 ± 0.3 


1.92 


Extended 


A = 


18 





5.3 ±0.4 


6.22 





6.1 ±0.3 


2.95 





5.1 ± 0.4 


6.70 


SNLS 


B = 


8 


5.6 ±0.8 





9.72 


6.1 ± 0.9 





9.74 


5.8 ±0.8 





10.09 


Extended 


B = 


18 


3.8 ±0.4 





9.42 


4.0 ±0.4 





9.56 


3.8 ±0.4 





9.66 














*i 


*2 




*i 


*2 




Extended 


2-bin d 


17 


90 ±5.2 


1.2 ±0.10 0.64 


120 ±7.8 


0.88 ±0.14 


0.79 


86 ±5.0 


1.3 ±0.10 


0.67 



a The extended sample refers to the SNLS sample plus the external data described in § [6j 
b Units of xl0 — 14 SNeyr -1 Mg , where Mq refers to the current stellar mass, M ste ii a r- 
c Units of xl0- 4 SNeyr- 1 (M yr- 1 )- 1 . 

d Units of xl0~ 14 SNeyr -1 Afg , where Mq refers to the total formed stellar mass, M*. 
e A discrete DTD, equal to *i at ( < 420Myr, and *2 otherwise. 
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As our default SFH model, we choose the[Lj (j2008l ) update to the lHopkins Beacoml (1200a) fi t 
to a compilation of recent star-formation density measures. For simplicity we use a lSalpeterl (|1955l ) 
initial mass function (IMF) with mass cut-offs at O.IMq and IOOMq, and assume that stars began 
to form at z = 10. The default SFH i s parameteriz e d in a piec e- wise fashio n , and is over-plotted 
on the data in Fig. [TTJ As shown by iForster et al.l (|2006l ) and iGraur et al.l (|201ll ) , the choice of 
SFH can add an additional significant systematic uncertainty in any comparisons of DTDs to SN la 
volumetric data. We the refore co mpare with results obtained using an alternative parameterization 
of the SFH byE il (l2008h followin g ICole et al.l (hoOlh. as well as a SFH fit to slightly different data 
by lYiiksel et al.l (|2008l ) (see also iHoriuchi &: Beacoml 1201 Oi l . In § 16.31 we also investigate a SFH 



derived in a very different manner (jWilkins et al 



200. 



The integral of the DTD gives the total number of SNe la per formed stellar mass, Ni a /M*. 
This can be converted into the fraction of intermediate-mass stars that explode as SNe la, 77, by 
multiplying by a factor of 47.2 (for the Salpeter IMF). This assumes that the progenitor mass range 
for a SN la is 3-8M Q (see iMaozl I2008I . for a discussion). For all our model DTDs, we set the DTD 
to zero at epochs earlier than 40Myr, the approximate lifetime of an 8M star. 



6.2.1. Gaussian DTDs 



We begin by fitting a Gau ssian DT D, with \fr(t) oc e~ ( - t ~ T ^ 2 ^ 2a2 \ to the volumetric SNRi a data, 
following IStrolger et al.l (|2004h (see also lStrolger et al.ll2010l ). We fit a DTD with parameters fixed 
at r = 3.4 Gyr and a = 0.2r (i.e., just adjusting the normalization in the fits), as well as a DTD 
fit with t allowed to vary. The results are listed in Table [4] and compared to other DTD fits in 
Fig. M 

This model has xi = 2.62 (xi 1S the reduced x 2 > the x 2 P er degree of freedom, v). Allowing r 
as a free parameter in the fits gives r = 3.1 ± 0.3 Gyr with a similar xl> ~ the fit quality is slightly 
better when using Cole et al. form of the SFH (x 2 , = 1.60). 



The Gaussian DTDs therefore provide poor fits to the SNRi a data, and are not capable of 
matching the SNLS, SDF and SDSS/LOSS data simultaneously. In particular, these Gaussian 
DTDs predict a decrease in the number of SNe at z > 1 not seen in the combined data set. However, 
they were originally favor ed following fits to data including z > 1 points from HST searches 



(Dahlen et al. 



2004 



20081 ) which are not included in our analysis due to their lower statistical 



pre cision compar e d to t he SDF study. As a consistency check we also replace the SDF data with 
the iDahlen et al.l (|2008l ) data (adjusted downwards by 15% to account for 19bg-lie events) in our 
fits - we find that the xi does not improve (Table H|) . 
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Fig. 13. — SN la rates as a function of redshift with various delay-time distribution (DTD) model 
predictions fit to the data for different cosmic SFHs. L ower left: U A + B" m odel (§ I6.2.3p . lower 
right: Gaussian DTDs (§ !6.2.ip . upper left: the model of lPritchet et al.l (|2008l ) (§ 16.2. 2p . and upper 
right: a gen eric p ower law DTD (§ 16.2.21). In a ll case s, solid lines repre sent t he piece- wise cosmic 
SFH of El (booA. dashed lin es the ICole et all ifcoOlh form of the <|200ah SFH, and dot-dash 
lines the IWilkins et al.l (|2008l ) SFH. See text for more details of the models, and Table H] for the 
numerical values. 
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6.2.2. Power law and exponential DTDs 



Theoretically, if SNe la are dominated by a single channel, the DTD will likely decline with 
age. In the single degenerate channel, SNe la at 10 Gyr should be rare, since 1 sec ondaries 



have small envelopes to donate and must rely on only the most massive primaries (jGreggid 120051 ) . 
A power law DTD (i.e., *$?(t) oc t~P with ~ 1) with a low time-delay cut-off is e xpected in the 
double degenerate scenario (e.g. lGreggidl2005l : lF5rster et al.ll2006l : iMaoz et alJl201Q ) . and has been 



expla ined post-hoc in the single degenerate channel using a mixture of contributions (jHachisu et al 
2008). Furthermor e, models with ~ 1 seem to provide a good match to a variety of recent 
observational data (jTotani et al.1 120081 : IMaoz et al.ll20ld ; iGraur et al.ll201ll ). 



We fit both = 1 and free DTDs to the SNLS+external SNRi a data. The results can 
be found in Table H and Fig. Q31 Our best-fit value is = 0.98 ± 0.05 (xl = °- 76 )> consistent 
with 1. This broad agreement with 1 holds when considering the other SFH paramterizations 
09 = 1.15 ± 0.08 for the Cole et al. form). 



Pritchet et al.l ([20081 . hereafter PHS) present a simple model relating white dwarf formation 



rate, which decreases with time following an instantaneous burst of star formation as ~ i~ 0,5 , 
resulting in a DTD with ~ 0.5. By fitting the SN la host galaxy data of IS06l . |PHS demonstrate 
that ^(t) ~ £-o.5±o.2^ i rreS p ec tive of the assumed SFH or the detailed mixture of stellar populations. 
PHS argue that the single-de generate forma tion s cenari o alone is not sufficient to account for all of 
the observed SNe la (see also lGreggidl2005l ). The iPHSl model makes an e xplici t prediction for the 
evolution of the SN la rate with redshift, given an input SFH. We fit the IPHSl model - essentially 
= 0.5 - to the data and show the resultant fit in Fig. [13] (also Table [3]). We find a xt = 4.08, 
obviously a substantially poorer fit than the generic power-law fit, or power-law DTDs with = 1. 

For completeness we also also test an exponential DTD, i.e., ^f(t) oc exp - */ T . When r is small, 
this approximates a simple star-formation dependence, and when large, it approximates a constant 
DTD. The results are in Table EJ generally, single exponential DTDs provide poor fits to the data, 
but do still have acceptable x 2 . 



6.2.3. "Two- component" models 

Finally, we examine various two-component DTD models. The first is the popular "A + £?" 
model, a simple, two-component model of SN la production that is comprised of a "prompt" 
component that tracks the instantaneous SFR, and a "delayed" (or "tardy ") component that is 
proportional to M ste u ar (jMannucci et al.ll2005l ; IScannapieco &; Bildstenll2005l ): 



SNR Ia (z) = A x M stellar (z) + B x SFR(z) 



(8) 



Here, the A and B coefficients scale the M ste n ar and SFR components, respectively. The prompt 
component consists of very young SNe la that explode relatively soon (in the model, immediately) 
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after the formation of their progenitors, whereas the delayed component (scaled by A) corresponds 
to longer delay times and an underlying old stellar population. This model is empirically attractive 
due to the ease of comparison with readily observable galaxy quantities, such as M s t e iiar and SFR. 
Note that this A+B model does not exactly correspond to a DTD - however, it can be easily 
converted to a DTD by using the variation of M ste n ar with time in a SSP. This leads to a DTD with 
some fraction of SNe la formed immediately (the B component), followed by a slightly decreasing 
fraction to large times (the A component). This decrease is ~ 25% from 0.1 to 5Gyr, and ~ 20% 
from 1 Gyr to 10 Gyr - clearly significantly shallower than a /3 = 1 power law. 

Some confusion exis ts over the exact definition of M ste u ar in eqn. ([8|), and hence the definition 



of A. Some authors (e.g. lNeill et al.ll2006l ) simply treat M ste u ar as the integral of the SFH, equating 
it to the total formed mass, M*. Others make corrections for stars that have died, particularly 
in studies which perform analyses on a gal axy-by-galaxy basis a s this quantity is more straight 



forward to link to observational data (e.g., ISullivan et al.l l2006bh . This latter definition leads to 



larger A values, as M s t e ii ar (i) will be less than M*(t) (see Fig. 7 of lSOa . for the size of this difference). 



Here, our A values ref er to M ste n ar , and we pass the cosmic SFH through the PEGASE.2 routine 



(|Le Borgne et al.ll2004l ). convolving the chosen SFH with a single stellar population, and generating 
a galaxy SED(i) from which mass M ste u ar (i) can be estimated. The evolving M ste u ar and SFR are 
used to perform a fit of equation [8] to the volumetric rate evolution, with results listed in Table |H 

Fitting the SNLS rates alone gives coefficients of A = (1.6 ± 0.5) x 10" 14 SNeyr" 1 Mq 1 and 



B = (3.4 ± 0.3) x lO^SNeyr™ 1 (A^yr^ 1 )- 1 for the QJ l|200d ) piece-wise SFH, with xl = 0.74 
for v = 7. Incorporating the external SN la rate from LOSS, SDSS and SDF yields values of 
A = (1.9 ±0.1) x lO-^SNeyr^M" 1 and B = (3.3 ± 0.2) x 10" 4 SNeyr" 1 (M yr" 1 )" 1 , with 
xl = 0.60. This fit is compared to other DTDs in Fig. [T3l 

Next, we set A = to investigate the possibility of a pure star-formation dependence, fitting 
only the prompt (B) component to the SNLS data. This results in an upper limit of B = (4.4 ± 
0.3) x 10 -4 SNeyr" 1 (Mq yr -1 ) -1 (xl = 1-81), equivalent to the normalized SFH curve plotted 
in Fig. [TTl Adding the external data rates to the SNLS values gives an upper limit of B = 
(5.3 ± 0.4) x 10~ 4 SNeyr -1 (M yr" 1 )" 1 , but again with a very poor fit quality (xl = 6.22). While 
the SNLS results themselves are marginally consistent with a pure prompt component (Table U]), 
adding the additional constraints supplied by the low-z data yield a very poor A = fit. The 
related test of setting B = and testing for only the delayed component similarly gives very poor 
fit results (Tabled!). 

As discussed above, these A and B values depend on the adopted SFH. Even ignoring the sys- 
tematics in the individual SFR measurements to which the SFH model is fit, the type of fit used also 
introduces considerable uncertainty. Using the Cole et al. form of the SFH, for example, changes the 
best-fit values to A = 1.5±0.2x 10~ 14 SNeyr^M" 1 and B = 4.3±0.3x 10" 4 SNeyr" 1 (il^yr^ 1 )- 1 , 
a significant variation with similar quality fit than for the piece-wise form. Thus care must be tak- 
ing in comparing any particular A + B values, or any particular prediction of SN la rate evolution, 
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without ensuring the consistent use of a SFH and derivation of M ste u ar . 



Clearly, and as well documented in the literature, the A + B model must be a significant 
approximation to the physical reality in SN la progenitor systems. In particular, there must be 
some delay time for the prompt component, and it is unlikely to act as a delta function in the 
DTD. We test this by approximating ^f(t) as two discrete bins in time, i.e., a step function with 
a value at tim es t < t sp n t , and ^2 at t > t S piit- We choose t sp m = 420 Myr (following, e.g., 
Brandt et al.ll2010l ). and used a sigmoid function to ensure the DTD was continuous when crossing 
itsplit- This DTD also provides a good fit to the data (Table 0]), with a significant detection of the 
two components - ^1 and '3/ 2 are > at 5<r in all three SFHs considered (and typically ~ 10cr). We 
experimented with making this function more general, by allowing t sp n t to vary. However, these 
fits were not constraining, although they prefer i S piit ^ 2 Gyr, and the fit parameters are highly 
correlated (i.e., as itsplit decreases, \&i increases to ensure a similar fraction of SNe la are generated 
from the "prompt" component). 

One direct outcome of this simple two-bin DTD model is that, for our default cosmic SFH, 
while ~ 70% of SNe la originate from the prompt component integrated over cosmic time, at z = 
the prompt component accounts for only ~ 25% of SNe la. These fractions remain fairly constant 
out to t S piit ~ 2 Gyr. 



We also exp lored more bins in \P - e.g., the three bin DTD of iBrandt et al.1 (l2010l ) and 
Maoz et al.l (|201ll ) - but, while these were consistent with our data, they did not provide improved 
fits over the two-bin DTD, and again, the parameters themselves were not well constrained. 



6.3. Discussion 



We now compare our DTDs inferred from the volumetric SNRi a data with other independent 
estimates from the literature, comparing both the parametric form of our best-fit DTDs, as well as 
the normalization. In Fig. [T4"l we plot our inferred DTDs, with normalizations from the best-fits to 
the volumetric SN la rates, and com pare to other empi r ical determinations of t he DTD from the 
literature from a variety of methods ( Brandt et al. 2010 : Maoz et al. 201ol 12011 ). The data points 
come from the analysis of SN la rates in galaxy clusters ( Maoz et al. 2010). the re construction of 



the DTD from an analysis of the LOSS SN la host galaxy spectra (jMaoz et al.ll201ll ). and a similar 



analysis of the host galaxies of SDSS SNe la (IBrandt et al. 



2010|). 



Where n ecessary we cony ert these external measure ments to a Salpeter IMF 
Salpet er IMF (jBell et alJboO^ fori 



from a diet- 



(|2010h . We also adjust the 



Maoz et al 



Maoz et al 



( 2010h and Maoz et al 



2011 



(|2010 . 2011) and a lKroupal (120071 ) IMF for lBrandt et al 



results downwards as their 
measurements presumably include SN1991bg-like events which are not included in our SNLS anal- 
ysis (the Brandt et al. analysis does not include these SNe, and t heir stretch range is well matched 
to our analysis - see their figure 2). Furthermore, we correct the IBrandt et al.l (|2010l ) DTD points 
upwards by 0.26 dex to account for stellar masses in Brandt et al. that are 0.26 dex too high due to 
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a normalization issue in the VESPA code used to derive them (D. Maoz, private communication; 
R. Tojeiro, private communication). 

Although the generic shape of the power-law DTD inferred from the volumetric rate data 
matches the external DTD data well (Fig. fT4"l) . it is clear that the best-fit normalization required 
to reproduce the volumetric rate data differs to that required to fit some of t he external samples . 
The DTDs inferred from volumetric rate da ta are genera l ly con sistent with the lBrandt et al.l ((2010) 
analysis, and the first and third bins of the lMaoz et al.l (|201ll ) data. 



However, the normalization of the best-fit DTD to the iMaoz et al.1 (|2010l ) cluster data lies 
significantly above our best-fit DTD. Integrating our best-fit power-law DTDs gives Ni^/M* ~ 
4.4 ± 0.2 - 5.2 ± 0.2 x lO^SNeM^ 1 (n = 2.0 - 2.5%) depending on the SFH, in good agreement 
with similar analyses (jHoriuchi &: Beacoml l2010l ). However, thi s is s ignificantly below the value 
of ~ 40 x 10~ 4 SNeMQ 1 obtained by integrating the lMaoz et al.l ( 20101 ) "optimal iron-constraint" 



power-law DTD (for our IMF), or ~ 24 x 10 -4 SNeM" 1 from the "minimal iron constraint" DTD. 
We can sanity check our normalizations by predicting the SN la rate in the Milky Way, given our 
DTD values. Assuming a Milky Way stellar mass of ~ 5 x 10 10 Mq and a SFR of ~ 4M Q yr~ 1 , 
our A + B DTDs give a predicted Milky Way normal SN la rate of 0.22-0.25 events per century. 
This is in go od agreement with independent estimates of the actu al rate (~ 0.35 to 0.40 events per 
century; e.g.. lTammannlll994l : iRoland et al.ll2006l : iLi et al.ll2011al ) given the uncertainties involved. 



So why is the normalization in the IMaoz et al.l (|2010l ) cluster SN la DTD (which is not arbi- 
trary) different to those derived from volumetric rate data by such a large factor? Several possibil- 
ities exist that may explain the discrepancy. The first is that the SNLS, and by extension all other 
SNRi a studies, are missing a significant number of SNe la, a factor of at least four. H owever, it is 



difficu lt to understand how this might occur, given some of the cluster rates used in IMaoz et al 
(J201CJ) are drawn from very similar surveys (including SNLS and SDSS-SN). 



A second possibility is that the cosmic SFH models used in our analysis over-predict the 
actual SFR - a lower SFH normalization would require a higher DTD normalization to match the 
volumetric rate data . To test this possibility, we repeat our rates fitting analysis using the SFH of 
Wilkins et al.l (|2008l ). a SFH derived from a requirement to match the redshift evolution of M ste iiar- 
Although this agrees with other SFH estimates at z < 0.7, it suggests that high-redshift SFRs could 
be ~ 0.6 dex lower compared to the models used in this paper. However, even using this SFH, the 
integrated power-law DTD only gives 6.6 ± 0.6 x 10~ 4 SNe Mg 1 (rj = 3.1 ± 0.3%), still some distance 
short of that apparently required by the clusters analysis. 

Other options to adjust the SFH normalization downward, such as reducing the dust extinc- 
tion corrections applied to the various star-formation indicators which make up the SFH compila- 
tions, are probably not viable given the long-established and significant evidence for obscuration 
in star- forming galaxies, an d the agreement between the different diagnostics (see discussion in 



Hopkins Beacom 



2006 



Lj|200a). 



A third possibility is that the cluster rates used to derive the DTD of lMaoz et al.l (|201dl ) have 
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Fig. 14. — The SN la delay-time distribution (DTD) inferred from fits to our volumetric rate data, 
compared to other determinations from the literature. The solid line shows the best-fitting power- 
law (PL) DTD and the long-dashed line the best-fitting "A + B" (AB) model from Table Efl each 
drawn for the two SF Hs that give the m ost different results. The short-dashed line is the best-fit 
power-law DTD from iMaoz et all (201C), which has /3 = 1.3, but with their normalization. The 
A + B model has been adjusted, for plotting purposes, so that the instantaneous c omponent is 
spread over the fi rst 400 Myr. The D TD data po i nts co me from iBrandt et al.l (|2010l . red circles), 



Maoz et al.l (|2010l . open circles), and IMaoz et al.l (|201ll . black circles). The horizontal error-bars 



indicate the bin widths on these points. 
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some contamination from "younger" SNe la, thus increasing their rates above that appropriate for 
their age, assuming a redshift of formation of 3 . This, and other s imilar potential systematics from 
the clusters analysis, are discussed in detail in lMaoz et al.l (|20ld ). 



Finally, it may well be the case that the assumption of a single DTD is not adequate, given 
the various indications that there may be more than one progenitor channel. This would suggest 
that there is not one universal DTD that is independent of redshift, or other variables, such as 
metallicity. 



7. Stretch dependence 

There is a n observed correlat i on between the photometri c properties of SNe la and their host 
environments (jHamuv et al.l Il995l . [200CJ; I Sullivan et al.l l2006bl ) . Brighter SNe la with slower light- 
curves tend to originate in late-type spiral galaxies, such that the rate s of higher-stretch objects are 
proportional to star- formation on short (~ 0.5 Gyr) timescales (jSOff). Meanwhile, fainter, faster- 
declining SNe la are more likely to be associated with older stellar populations. This split seems 
to extend to the recovered DTDs - iBrandt et al.l (j2010l ) show that the recovered DTD for low and 
high stretch SNe la are very different, consistent with the above picture of young SNe la being 
high stretch and old SNe la low stretch. A larger fraction of high-stretch SNe la might therefore 
be expected at high re dshift tracking the increase in the cosmic SFH and hence preponderance of 
younger stars. Indeed, iHowell et al.l (|2007l ) find a modest increase in the average light-curve width 
out to redshifts of z ~ 1.25. The rate evolution for low-stretch SNe la should therefore demonstrate 
a correspondingly shallower increase with redshift than that of higher-stretch objects. 

We investigate this trend in Fig. [T5l splitting the SNLS sample at the median stretch value of 
so = 1.0. The measured rates with simple weighted errors are plotted separately for objects with 
0.8 < s < 1.0 (open cir cles) and those with 1.0 < s < 1.3 (filled circles). We also show the s < 0.8 



rate measurement from lGonzalez-Gaitan et al.l (|201ll ). The samples from the analysis in this paper 
exhibit a comparable rise in their rates with redshift, with power-law slopes of ao.8<s<i = 2.48±0.97 
and a s >i = 2.11 ± 0.60. Extrapolating the fits to each of the samples back to z = gives the 
following fractions of SNe la at z = in each group: 17% {s < 0.8), 39% (0.8 < s < 1.0), 



fraction of 15% ( 


Li et al. 


2011a 


( Gonzalez-Gaitan et al. 


201 1|). r 



2011). The fraction of very lo w stretch events with s < 0.8 SNe la shows 



uncertainties are very large. 

The stretch-split rates are only considered out to z = 0.8, beyond which redshift the stretch 
errors become large (> 0.1), and the lower efficiencies of the redder, lower-stretch SNe la can bias 
the observed stretch evolution by driving up the s < 1 rates (Fig. [7|). 



Fig. [HI shows the ratio of the rates split by stretch for s > 0.8. To compare this observed 
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Fig. 15. — SN la rates split at the median SNLS stretch value of so = 1-0. The lower-stretch rates 
are shown as open circles, and high-stretch rates by the solid circles. We additionally show the 
s < 0.8 rate measurement from iGonzalez-Gaitan et al.l (|201ll ). The measured rates have simple 
weighted errors and are uncorrected for systematic errors in redshift. The first two redshift bins 
{z = 0.15,0.25) are combined here due to the small number of low- z objects the split samples. 
Power-law fits to the rates yield similar slopes over the redshift range shown: a s <o.s = 1-63 ± 0.28 
(dotted line), a .8<s<i = 2.48 ± 0.97 (dashed line) and a s >i = 2.11 ± 0.60 (solid line). 
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Fig. 16. — The ratio of low-stretch (0.8 < s < 1) to high-stretch (1.0 < s < 1.3) SNe la volumetric 
rates as a function of redshift. The horizontal line shows the weighted mean ratio of all SN la 
stretc hes in the sample of 1.01. Various predicted trends, following the analysis of iHowell et al 



(|2007l ) and using the A + B and power-law DTDs, are overlaid. These trends are not fit to the 
data-points plotted. The solid lin es rep resent the piece-wise SFH, and the dashed lines the Cole 
et al. form of the SFH, from Q ( 2008 ). The shaded areas show, for each pred icted trend, the 
uncertainty expected by shifting the stretch distributions of IHowell et al.l (|2007l ) by 0.025. The 
reduced % 2 of the model fits to the data range from xt = 1-8 to xt = 2.6 over the redshift range 
shown, compared with x„ = 1-35 for a flat line at the weighted mean ratio. 
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data with any expected evolution, we need model stretch distributions for "old" and "young" SNe, 
together with a mechanism for predicting the relative evolution of these two components with 
redshift. For the fo rmer, we take the two stretch distributions for young and old SNe la from 



Howell et al.l (|2007l ): the old component is represented by a Gaussian with (s ld) = 0.945 and 
cr Sold = 0.077, while the young component has (s yoU ng) = 1-071 and cr Syouns = 0.063. To estimate 
the relative redshift evolution of the old and young components, we use the best-fitting A+B values 
from Tabled] (assigning A to the old SNe and B to the young SNe), as well as the power-law DTD. 
For this latter DTD, we assign SNe born at t < 2Gyr in the DTD to the young component, and 
SNe born at t > 2Gyr to the old component. Together with a SFH, these models then predict 
the relative fracti on of low and high stretch SNe as a function of redshift, over-plotted in Fig. [T6j 



We also vary the lHowell et al.l (|2007l ) stretch distributions by adjusting (s) by ±0.025 for the two 



components; these are shown as the hashed gray areas in the figure. 

As expected, the predicted ratios show a smooth decline from a large fraction of low-stretch SNe 
at z = 0.1. As the relative contribution from delayed SNe la to the rates decre ases with increasin g: 
redshift, so too should the dominance of lower-stretch objects (see also Fig. 1 in lHowell et al.ll2007l ). 
The prediction based on the power-law DTD shows a shallower evolution with redshift, reflecting 
the extended age of the young component relative to the simplistic A + B model. 

However, Fig. [15] and Fig. [16] are surprising. If broad- lightcurve SNe la favor a young environ- 
ment and narrow-lightcurve events favor an old environment, as has been well established, then the 
ratio of narrow to broad SNe ought to be changing as star formation increases with redshift. But 
all of the predictions are a relatively poor match to the SNLS sample in the higher redshift bins. 
We vary so by ±5% to assess the sensitivity of our results to the stretch split value (default of 1.0), 
but find no significant improvement in the agreement with the predicted model as compared with 
the straight-line fit at the weighted mean ratio. 

The lack of an observed evolution may be due to several factors. The first is that the age-split 
between low- and high-s SNe la may be more subtle than previously appreciated. Alternatively, 
the main effect may be dominated by s < 0.8 SNe la, which do show a different evolutionary trend 
with redshift (Fig. 1150 . The limited time baseline of only ~ 4.5 Gyr from z = 0.2 — 0.8 may also be 
a factor, and of course limitations of the method (e.g., the arbitrary cutoff of for SNe to be "young" 
or "old" in the power-law DTD) could mask any real effect. 

Some other aspects of Fig. [16] are better understood. The A + B model over-predicts low-s 
SNe at z = because it has only 20% fewer SNe in the DTD at 12 Gyr compared to 1 Gyr, in 
apparent contrast to the data and to the t _1 model, where the rate falls by an order of magnitude 
over this baseline. These excess old SNe la from z = 1.5 show up after a 10 Gyr delay at z = 0. 
However, the difference in the predictions between the A + B and power-law models is not large, 
so this unphysical assumption cannot provide the only solution for the lack of observed evolution 
in Fig. [16] 
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8. Summary 



In this paper, we have probed the volumetric rate evolution of "normal" 0.8 < s < 1.3 SNe la 
using a sample of 691 events from the Supernova Legacy Survey (SNLS) in the range 0.1 < z < 
1.1, 286 of which have been confirmed spectroscopically. The SNLS rates increase with redshift 
as (l+z) a with a = 2.11 ± 0.28, and show no evidence of flattening beyond z ~ 0.5. Due to 
spectroscopic incompleteness and the decrease in detection efficiency for the SNLS sample, a rollover 
in the slope cannot be ruled out beyond z ~ 1 based on the SNLS data alone. 

As a significant component of the SN la rate is linked with young stellar populations, an 
increasing fraction of SN la events may suffer the effects of host extinction at higher redshifts. In 
our rate calculation method, the effect of SN color is factored directly into the detection efficiency 
determinations: detection recovery is evaluated empirically according to the observed SN color 
regardless of its cause. Redder objects at a given redshift have lower detection efficiencies, and are 
correspondingly more heavily weighted in the rates determination. 

Combining the SNLS data with that from other SN la surveys, we fit various simple delay- 
time distributions (DTDs) to the volumetric SN la rate data. DTDs with a single Gaussian are not 
favored by the data. We find that simple power-law DTDs oc t _/3 ) with /3 ~ 1 ((3 = 0.98±0.05 

to (3 = 1.15±0.08 depending on the parameterization of the cosmic SFH) can adequately explain all 
the SN la volumetric rate data, as can two-component models with a prompt and delayed channel. 
These models cannot be separated with the current volumetric rate data. Integrating these different 
DTDs gives the total number of SNe la per solar mass formed (excluding sub-luminous s < 0.8 
events) of N\ a /M* ~ 4.4 — 5.7 x 10~ 4 SNeMg 1 (assuming a Salpeter IMF), depending on the star 
formation history and DTD model. This is in good agreement with other similar analyses, but lies 
significantly below the number expected from DTDs derived from cluster SN la rates. 



The use of other techniques, such as fitting the SFH of individual galaxies (jSullivan et al.ll20 06b; 



Brandt et al 



200 



2010 



Maoz et al.1 



Maozeta! 2011 



(Sullivan et al. 


2006b: 


galaxies ( 


Totani et al. 



20 id ), use more information, and in principle oug ht to be more reliab l e. However, 



each technique has significant drawba cks, such as contamination ( Totan i et al 



20111 ). limitations of SED fitting codes (ISullivan et al.ll2006bl : iBrandt et al.ll20l"n[lMaoz et al.ll201ll ). 
and the assumptio n that all cluster galaxies formed at z = 3 in a delta-function of star-formation 



2008 



Maoz et al 



(|Maoz et al.ll201ll ). Therefore, our results are an important complementary constraint. By pre- 
senting an evolution in the SN la rate over a large redshift baseline done self-consistently by a 
single survey we have for the first time mitigated the primary drawback of this method - having 
to combine myriad rate determinations from mul tiple surveys, all done with different assumptions 
and biases, sometimes disparate by large factors (jNeill et al.l l2006). 



We also find no clear evidence for a difference in the rate evolution for SNLS samples with 
0.8 < s < 1.0 and 1.0 < s < 1.3 out to z = 0.8, although the stretch evolution model from 



Howell et al.l ([20070 cannot be ruled o ut conclusively. Stretch evol ution plays a more significant 



role in the sub-luminous population ([Gonzalez-Gaitan et al.l l201ll ) , which show a much flatter 
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evolution than the s > 0.8 sample. 

Next generation surveys such as Dark Energy Survey (DES), Pan-STARRS, Palomar Transient 
Factory (PTF), and SkyMapper, many of which are already underway, are finding thousands of SNe 
la (in comparison to the ~ 700 in this study). Statistical rate determinations ought to improve, 
but systematic difficulties will remain, as not all SNe can be spectroscopically confirmed. However, 
large number statistics will allow the construction of sub-samples larger than the three (split by 
stretch) analyzed here. Comparison of the relative rates of SNe with different properties and in 
different environments may ultimately improve deduced DTDs, and allow for the construction of 
different DTDs for subsets of SNe la. 
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